Orbital ordering in the ferromagnetic insulator Cs 2 AgF 4 from first principles 
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We found, using density-functional theory calculations within the generalized gradient approxi- 
mation, that Cs2AgF4 is stabilized in the insulating orthorhombic phase rather than in the metallic 
tetragonal phase. The lattice distortion present in the orthorhombic phase corresponds to the 
x 2 — z 2 /y 2 — z 2 hole-orbital ordering of the Ag 2+ 4d 9 ions, and this orbital ordering leads to the 
observed ferromagnetism, as confirmed by the present total-energy calculations. This picture holds 
in the presence of moderate 4d-electron correlation. The results are compared with the picture of 
ferromagnetism based on the metallic tetragonal phase. 

PACS numbers: 71.20.-b, 75.25.+Z 



I. INTRODUCTION 

It is quite often that the 3c? transition-metal ox- 
ides with orbital degeneracy display an orbital order- 
ing associated with lattice distortions (of the Jahn- Teller 
type). 1 ' 2 In contrast, -id and 5d transition-metal com- 
pounds rarely show this^^^i 7 - mainly due to the de- 
localized character of the d electrons. It was therefore 
a surprise when the 4c?-orbital ordering and the orbitally 
driven spin-singlet dimerization were found in the ruthen- 
ate La4Ru20io^ Moreover, the fascinating superstruc- 
ture of the spinel CuIr2S4 (Ref. 0) was explained in terms 
of the concept of orbitally-induced Peierls state (associ- 
ated with the partially occupied Ir 5c? orbitals)^ 

Very recently;^ the magnetic behavior of the layered 
Cs2AgF4 was investigated, and the in-plane ferromag- 
netism (FM) of the 5=1/2 Ag 2+ ions was evidenced by 
magnetization and inelastic neutron scattering measure- 
ments, in sharp contrast to the structurally analogous 
high-T c cuprates. Moreover, analysis of the structural 
data suggests that Cs2AgF 4 is 4c? orbitally-ordered, and 
it was proposed that such orbital ordering may be the 
origin of the observed FM. 

Cs2AgF4 was first synthesized in 1974 by Odenthal, 
Paus, and Hoppe^ and it received, as well as other 
silver fluorides, attention due to pursuit of supercon- 
ductivity in transition-metal compounds other than the 



cuprates 
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Cs2AgF4 is a structural analog of the 



parent compound of cuprate superconductors, La2Cu04 
(AgF2 vs CuC>2 sheet), and both have the 5=1/2 diva- 
lent transition-metal cations (Ag 2+ 4c? 9 vs Cu 2+ 3c? 9 ) and 
isoelectronic anions (F~ vs 2 ~). It is well known that 
La2CuC>4 is a Mott insulator of charge-transfer type with 
two-dimensional Heisenberg antifcrromagnetism (AF). 
Upon hole doping (like Ba substitution for La) it is tuned 
into a superconductor, and, its high-temperature normal 
state being a paramagnetic metal, has the holes predom- 
inantly on the antibonding pda band. 14 Silver fluorides 
have similar electronic band structure, namely the Ag-F 
pda antibonding character in the vicinity of the Fermi 
levels Moreover, the Ag-F bonding is substantially co- 
valent in the Ag 2+ and Ag 3+ fluorides, and thus holes 
might be doped into the F 2p band, which is unprece- 



dented in the chemistry of transition-metal compounds. 
As a result, an interesting new perspective is opened for 
the design of a novel family of metallic, perhaps even 
superconducting materials based on silver fluoridesji 2 - 

Subsequent to the finding of FM in Cs2AgF 4f £ its elec- 
tronic structure was studied by Kasinathan et al^ and 
Dai et al.^ both using density-functional theory (DFT). 
Kasinathan et al. found the FM, due to importance of the 
Ag-F covalency, to be itinerant in character, and substan- 
tial magnetic moment on the fluorine ions as result of a 
Stoner instability enhanced by Hund's coupling^ Dai et 
al. suggested that the FM originates from the spin po- 



larization of the doubly occupied 



y band, which 



is induced by the c? 2 2 — p — d x i_ y i orbital interaction 
through the Ag-F-Ag bridges.^ Both groups of authors 
dealt with the tetragonal lattice 10 and obtained a nearly 
half-metallic solution in the FM state, and used it to 
explain the origin of FM in Cs2AgF4 i 15 ' 16 However, in 
reality Cs2AgF 4 is an insulator with a relatively large 
energy gap (lilac color) and has a distorted orthorhom- 
bic lattice,— and thus the above studies of the FM are 
insufficient. More recently, Kan et alJ^ performed pseu- 
dopotential DFT calculations for the orthorhombic struc- 
ture and they obtained an orbitally ordered solution (in 
contrast to Refs. [l5l and [Tfjh . However, their calculations 
underestimated the in-plane distortion and strongly over- 
estimated both the in-plane and out-of-plane magnetic 
coupling strengths. Moreover, they did not investigate 
the tetragonal structure. Therefore, their results do not 
allow to reach solid conclusions. 

In the present work, we re-investigate theoretically the 
electronic structure of Cs2AgF4 and the origin of its FM, 
by carrying out a comparative study of the orthorhom- 
bic and tetragonal lattices of Cs2AgF4 using DFT total- 
energy calculations. We find that the orthorhombic phase 
is more stable and it is insulating with suppressed Ag-F 
bandwidths. The inherent lattice distortion is accom- 
panied by the Ag 4c?-orbital ordering, and this orbital 
ordering readily accounts for the observed FM (similar 
results were also obtained recently by Hao et alJ£). This 
picture holds in the presence of a moderate correlation 
of the Ag 4c? electrons. The resulting picture of an or- 
bital ordering and of a ferromagnetic insulating state of 
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Cs2AgF4 is thus very similar to that of an isoelectronic 
and isostructural compound K^CuF^^ 



II. RESULTS AND DISCUSSION 

Our electronic structure calculations are performed 
by using the all-electron full-potential augmented plane 
wave plus local orbital method^ The generalized gradi- 
ent approximation^! (GGA) to density-functional theory 
is adopted.— We took the neutron diffraction structure 
data measured at 6 As seen below, the electronic 
structure and magnetic properties of Cs2AgF4 are de- 
termined by an orbital ordering, and the orbital order- 
ing is the most relevant to the local distortion of the 
AgFg octahedra. To account for a possible lattice dis- 
tortion connected with an orbital ordering, we carried 
out an optimization of the atomic positions, keeping the 
unit-cell parameters fixed and relaxing the atomic coor- 
dinates. The obtained values of atomic positions and 
bond-lengths are in good agreement with the experimen- 
tally determined values, see below. Based on that, in 
the further calculations we use the experimental struc- 
ture data. The Cs hphdQs (4d5s), Ag Ad5s5p (4s4p), and 
F 2p3s (2s) are treated as valence (semicore) states. The 
muffin-tin sphere radii are chosen to be 2.8, 2.2, and 1.7 
Bohr for Cs, Ag, and F atoms, respectively. The cut- 
off energy of 16 Ryd is set for the plane-wave expansion 
of the interstitial wave functions. The sufficiently dense 
k mesh, e.g., 10x10x4 is used for integration over the 
Brillouin zone of the primitive cell with the experimen- 
tal orthorhombic lattice constant a=6.4345 A , 6=6.4390 
A , and c=14.1495 A & The present set-up ensures a 
sufficient accuracy of the calculations. 

Fig. 1 shows the total and the orbital-resolved den- 
sity of states (DOS) for the FM orthorhombic phase of 
Cs2AgF4 obtained in GGA. It is an insulator with a gap 
of about 0.2 eV, but neither a metal (nor a half-metal) 
as predicted before for the tetragonal phas o 15 ' 16 (see also 
below) nor a semimetal as obtained for the orthorhombic 
phased The strong hybridization between Ag Ad and 
F 2p levels (apical F ap and planar F p i), both centered 
around 3 eV below the Fermi level, is obvious in both 
the valence and the conduction bands, and this gives rise 
to the large bonding-antibonding splitting of about 2 eV 
for the Ag t2 g (xz, yz, and xy) orbitals and even a larger 
one for the e g orbitals. The t% g bands are fully occu- 
pied and the e g bands are three-fourth occupied (with 
one hole left on them) as expected for the formal Ag 2+ 
(Ad 9 ) ions. Actually, the one hole spreads, due to the 
strong Ag-F covalency, over the six fluorine atoms of the 
AgF6 octahedron. As a result, the hole state consists of 
approximately 60% of the Ag Ad and 40% of the 2p of 
the four F atoms (see the DOSs above the Fermi level in 
the second, third, and fourth panels of Fig. 1), according 
to the stoichiometry of Cs2AgF4^ This is also reflected 
by the magnitude of local spin magnetic moments within 
each muffin-tin sphere, 0.560 /^s/Ag, 0.105 /^s/F ap , and 
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FIG. 1: (Color online) Density of states (DOS) of Cs 2 AgF 4 
in the FM orthorhombic phase from GGA. This insulating 
solution has a small band gap of 0.2 eV. The Fermi level is 
set at zero energy. Solid red (dashed blue) lines depict the 
spin-up (down) states. The panels show, from top to bottom, 
the total DOS per formula unit and orbital-resolved DOS. 
The state around -7 eV comes from the Cs 5p orbital. For 
the Ag 4<i-resolved DOS, the local coordinate system is chosen 
in such a way that the y (ir)-axis is along the in-plane longer 
(shorter) Ag-F bond. The Ag 2+ x 2 — z 2 hole state, and the 
corresponding y 2 — z 2 hole state of the nearest neighbor Ag 2+ 
ions form the in-plane orbital ordering (see also Fig. 2). 



0.112 jUs/Fpi, see also Table I. The hole states have al- 
ternating x 2 — z 2 and y 2 — z 2 symmetry shown in Figs. 
1 and 2. 

Thus, even the GGA calculations, not including corre- 
lation effects, already give an orbital ordering of the same 
type as that in K2CuF4^ which in a picture of local- 
ized electrons would immediately give in-plane ferromag- 
netism, according to Goodenough-Kanamori- Anderson 
rules. This orbitally-polarized hole state in the insulat- 
ing solution is associated with the large distortion of the 
Ag-F bonds of 0.33 A (2.44 A along the y-axis ver- 
sus 2.11 A along the a;-axis, and 2.11 A along the z- 
axis, all in a local coordinate system)^ which leads to a 
pronounced crystal-field splitting and suppressed band- 
widths. As a result, the antibonding x 2 — z 2 band lies 
about 1 eV above the antibonding 3y 2 —r 2 band, and then 
an exchange splitting of the x 2 — z 2 band gives rise to the 
small band gap and the x 2 — z 2 hole state. Thus, such a 
cooperative distortion accommodates the x 2 — z 2 jy 2 — z 2 
hole-orbital ordering in the AgF2 basal plane, giving the 
observed FM. Also, our GGA calculations confirm this 
by showing that the in-plane (intralayer) FM state is in- 
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TABLE I: Electronic structure of the in-plane FM and AF states of Cs2AgF4 in the orthorhombic and tetragonal phases 
calculated by GGA and GGA+U with U—3 eV. Energy difference (AE, meV/f.u.), gap size (eV), and local spin moments (/is) 
of each Ag, apical F (F ap ) and planar F (F p i) atom are shown. The AF solution of the tetragonal phase by GGA is close to a 
nonmagnetic (NM) solution indicated in the bracket. Note that the orthorhombic FM insulating solution is the ground state. 
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FIG. 2: (Color online) Contour plot (0.1-0.8 e/A 3 ) of the 
spin density in the the xz (upper panel) and yz (lower panel) 
planes, through the AgFo octahedra of Cs2AgF4 in the FM 
insulating orthorhombic phase obtained by GGA. It evidences 
the x 2 — z 2 /y 2 — z 2 hole-orbital ordering of the Ag 2+ ions in 
the xy basal plane. See also Fig. 1. Note also a relatively 
strong spin polarization of the F atoms (see more in the text). 
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FIG. 3: (Color online) Orbital-resolved density of states 
(DOS) of Cs 2 AgF 4 in the FM tetragonal phase from GGA. 
It is nearly half-metallic, and the wide in-plane x 2 — y 2 band 
and the small e 3 -level splitting are obvious. 



deed more stable than the in-plane AF state by 34 meV 
per formula unit (f.u.), see Table I. We note that the 
interlayer magnetic coupling is too weak (0.004 meV in 
energy scaled) to be captured by the present calculations, 
which show practically degenerate interlayer FM and AF 
states within the error bar of 1 meV/f.u. 

Actually, the experimental distorted lattice turns out 
to be in the vicinity of the theoretical equilibrium state 
after atomic relaxation. The structural optimization 
shows that Cs would shift by only 0.005 A and the Ag- 
F ap bond would be elongated by 0.041 A and the in-plane 
longer (shorter) Ag-F p ; bond would be shortened (elon- 
gated) by 0.018 A , as compared with the experimentally 
observed structure. As a result, the large in-plane dis- 
tortion of about 0.3 A , found experimentally,— is well re- 
produced theoretically. Moreover, the FM orthorhombic 
phase turns out to be more stable than a FM tetrago- 
nal phase by 27 meV/f.u. (which, as well as the follow- 



ing GGA±C/ result, implies a stability of the orthorhom- 
bic phase above room temperature as observed 8 ). We 
checked it by performing the calculations for the average 
tetragonal structure . 

The tetragonal phase is constructed in such a way that 
the lattice volume, c-axis constant, and positions of Cs 
and F ap atoms are the same as in the orthorhombic phase 
we calculated above, and that the basal plane is a regular 
AgF 2 square lattice. This FM tetragonal phase shows the 
same electronic-structure features (see, e.g., near half- 
metallicity in Fig. 3) as previously found ) 15 ' 16 since the 
structural parameters used in fiefs. [TH and [l6| are very 
close to what we used for this hypothetical phase. The 
regular AgF 2 square sheet would produce a wide x 2 — y 2 
band with a bandwidth of about 3 eV, and the com- 
pressed distortion of the AgFg octahedron (out-of-plane 
Ag-F bond length of 2.12 A vs in-plane one of 2.29 
A ) is not large enough to generate a crystal-field split- 
ting in excess of the bandwidths to open a gap. (Note 
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also that the compressed octahedron is extremely rare in 
Jahn- Teller insulators with e g degeneracy^) Actually, 
the crystal field splitting would make the x 2 — y 2 level lie 
below 2>z 2 — r 2 by only about 0.2 eV. As a result, both the 
minority-spin x 2 — y 2 and 3z 2 — r 2 bands cross the Fermi 
level. The total spin moment of 0.98 /is/f.u. (close to 
integer 1 /is ) , in the tetragonal phase would spread over 
Ag (0.558 fi B ), F ap (0.103 \i B each), and F pl (0.117 \i B 
each). 
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FIG. 4: (Color online) Density of states (DOS) of Cs 2 AgF 4 
in the FM orthorhombic phase from GGA+U with (7=3 eV. 
The insulating gap is increased up to 1.0 eV, compared with 
GGA (see Fig. 1). See other notes in Fig. 1 caption. 

As we saw above, Cs 2 AgF 4 is insulating in the FM- 
ordered orthorhombic phase, even in GGA, see Fig. 1. 
A moderate electron correlation may be present for the 
Ag 2+ Ad electrons, and thus we also carried out cal- 
culations using GGA+C7 method 25 to include the on- 
site Coulomb interaction for both the orthorhombic and 
tetragonal phases. In Fig. 4 we show the DOS results 
for the actual FM orthorhombic structure, given by the 
GGA+J7 calculation with an effective U—3 eV for the 
Ag 2+ Ad electrons. The insulating gap is increased up to 
1.0 eV (1.4 eV when £7=4 eV). The Hubbard U pushes 
the occupied Ad levels downwards so that their center of 
gravity is now lower in energy than the F 2p levels. As 
a result, the lower- lying bonding states have larger Ag 
Ad character than the higher-lying antibonding states, in 
contrast to the GGA results shown in Fig. 1. Again, 



the x 2 — z 2 /y 2 — z 2 hole-orbital ordering is obvious, and 
the orbital polarization of the e g states is enhanced. The 
local spin moments are 0.600 ^b/ Ag, 0.099 /is/F ap , and 
0.097 hb/F p i. Note that in the GGA+C7 with U=3 eV 
(U=A eV), the in-plane FM state is more stable than 
the in-plane AF state by 13 meV (10 meV) per formula 
unit. Using a simple Hcisenberg model for the spin-1/2 
square lattice, H = —JJ2(ij)SiSj (counting twice the 
magnetic exchange per spin pair), the nearest-neighbor 
intersite FM exchange integral is estimated to be 6.5 meV 
(5 meV) when U=3 eV (U=A eV). This value, reduced 
from the GGA estimate of 17 meV, is well comparable to 
the experimental one of 3.8-5 meVfS 

The GGA+U ([7=3-4 eV) calculations for the tetrago- 
nal structure still give a metallic solution for the in-plane 
FM state with the broad e g bands, but open an insulat- 
ing gap of 0.4-0.7 eV for the in-plane AF state with sup- 
pressed bandwidths (not shown). With the gap opening, 
this insulating AF state becomes stable with a large Ag 
spin moment, in contrast to the GGA results (see Table 
I). Note however that also in GGA+U both the AF and 
FM tetragonal phases have much higher energy than the 
orthorhombic phases, by more than 80 meV/f.u.. Thus, 
once again these calculations show that the orthorhombic 
phase is indeed more stable than the tetragonal phase, 
with the orbitally ordered FM state being the ground 
state. All this confirms the recent experimental finding 
that the layered Cs2AgF4 is indeed stabilized in the in- 
sulating orthorhombic phase. 8 



III. CONCLUSION 

Summarizing, by means of GGA and GGA+C7 band 
structure calculations, we find that the layered Cs2AgF 4 , 
a structural analog of cuprates, is stabilized in an in- 
sulating orthorhombic phase rather than in a metallic 
tetragonal phase. The intrinsic lattice distortion of the 
orthorhombic phase is accompanied by the x 2 — z 2 /y 2 — z 2 
hole orbital ordering, which readily accounts for the ob- 
served in-plane ferromagnetism. The present calculations 
confirm the recent experiments 8 - and lead us to the con- 
clusion that Cs2AgF4 is indeed an orbitally-ordered fer- 
romagnetic insulator, in strong contrast to the antiferro- 
magnetic nature of the parent high-T c cuprates, but in 
close analogy to K 2 CuF 4 . 
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